Casimir force between planes as a boundary finite size effect 
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The ground state energy of a boundary quantum field theory is derived in planar geometry in 
D+l dimensional spacetime. It provides a universal expression for the Casimir energy which exhibits 
its dependence on the boundary conditions via the reflection amplitudes of the low energy particle 
excitations. We demonstrate the easy and straightforward applicability of the general expression by 
analyzing the free scalar field with Robin boundary condition and by rederiving the most important 
results available in the literature for this geometry. 
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I. INTRODUCTION 

The Casimir force between two neutral macroscopic 
bodies in vacuum is considered as a manifestation of zero 
point fluctuations intrinsic to any quantum field theory 
(QFT). It has received a lot of interest recently both from 
the theoretical and from the experimental side (cf. the 
recent reviews Q and the book by Milton 0). 

This force can be determined from the volume depen- 
dence of the ground state energy (Casimir energy). It 
is usually derived either by a summation of zero point 
fluctuations of the quantum field restricted to the do- 
main in question, or by analyzing the statistical fluctu- 
ations of the interactions between the boundaries. Our 
new derivation is similar to the first approach in the sense 
that it uses the quantum fields living between the bound- 
aries. However, we describe the quantum field theory in 
question in terms of its particle excitations rather than 
the field fluctuations. As a consequence we describe the 
Casimir energy as 'virtual particle' exchange between the 
two boundaries and our formula is expressed in terms of 
the analytically continued reflection amplitudes of these 
particles. 

The interaction between the quantum field and the 
macroscopic bodies is most often described by imposing 
appropriate boundary conditions on the field. The aim 
of this paper is to explore and make explicit the depen- 
dence of the Casimir force on the boundary conditions 
obeyed by the field(s) in a geometry when the boundary 
conditions are given on two parallel planes. We achieve 
this by regarding the force as a physical manifestation of 
finite size effects in a boundary QFT (BQFT). 

In QFT, finite size effects in a periodic box were con- 
sidered by Luscher who showed how the correction 
to particle masses (infinite volume energies) can be ex- 
pressed in terms of infinite volume scattering data (5 1 
matrix elements) . Following this general idea, we have re- 
cently determined the finite size correction to the ground 
state energy in 2 dimensional BQFT in terms of scatter- 
ing data on the half line (reflection amplitudes). This 
explicitly describes how the ground state energy depends 
on the volume and it is obvious from the results in Q 
that a generalization to BQFT in any spacetime dimen- 



sion is straightforward and leads to a description of the 
Casimir effect. 

Boundary conditions in quantum field theory can be 
equally well described by a boundary state in a crossed 
channel. This formalism was worked out in detail for 
integrable 1+1 dimensional BQFTs |5|, but it can be 
extended to nonintegrable field theories in any number 
of spacetime dimensions. We show that this formalism 
allows an alternative method to derive the Casimir en- 
ergy, and it has the advantage of using only renormalized 
and phenomenologically meaningful field theoretic quan- 
tities. A further advantage of using the boundary state 
formalism is that it provides a systematic (large volume) 
expansion for the Casimir energy in the interacting case. 

Our main result is a large volume expression for the 
Casimir energy/force, valid in any BQFT, which depends 
only on the reflection amplitudes on the boundaries (al- 
though the boundary conditions considered allow trans- 
mission as well). These reflection amplitudes can be de- 
termined in a straightforward way in the half-infinite (one 
boundary) geometry, and using them in the expression 
for the Casimir energy makes it unnecessary to carry out 
the finite volume quantization separately in the various 
applications. 

The paper is organized as follows: in Section II we 
introduce BQFTs and derive our main results for the 
volume dependence of their ground state energy using 
the boundary state formulation. Section III contains a 
description of several applications of the result in vari- 
ous BQFTs: First we derive a formula for the Casimir 
energy of a free massive boson subject to Robin type 
boundary condition. This expression is tested against the 
well-known Dirichlet, Neumann and massless limits. We 
further check our formula by recomputing the Casimir 
effect for parallel dielectric slabs separated by a vacuum 
slot and for a massless fermion subject to "bag boundary 
condition" . We make our concluding remarks in Section 
IV. In the Appendix we confirm our result in a simple 
case by presenting the naive derivation based on mode 
summation tailor made for the framework of boundary 
QFTs. 
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II. DERIVATION OF THE MAIN FORMULA 

Here we summarize the main aspects of BQFTs and 
present the novel derivation of the Casimir energy. We 
follow the description of BQFTs along the line of 6] ex- 
tending the result of [f| for higher dimensions. In the 
process we show that this formalism makes it possible 
to establish a systematic (large volume) expansion of the 
Casimir energy even for interacting fields. 



A. Boundary state formalism 

To simplify the presentation consider a quantum field 
theory of a scalar field, $>(t,x,y), in D + 1 dimensions 
with Lagrangian 

c = \ - \ (a,*) 2 - ~ (aV) 2 - v@) 

where y is a D — 1 dimensional position vector, di — 
-p-. We also suppose that the spectrum consists of one 
particle type with mass m. The Hilbert space of the 
model can be spanned by asymptotic multi-particle states 
with momentum parameterized by (ki,ki) (where again 
ki denotes a D — 1 dimensional momentum vector) : 

H = fci; fc 2 , k 2 ; ■ ■ ■ ; k n , k n )\ 
The energy of a one-particle state is 



ui(k, k) = Y to 2 + k 2 + k 2 

In the asymptotic scattering configurations (large nega- 
tive or positive time) the particles are distant from each 
other, and the spectrum can be described by free in or 
out particles so the energy of the multi-particle state is 
the sum of the individual energies. These two Hilbert 
spaces are connected by the scattering matrix 

S({k i , fc^}, {ki, ki}) — 

(fcj, k-^', fej) &2) "•' ^n' 1^1 J ^1 i ^2 > &2 i k m , k m ) 

which can be expressed in terms of the correlators 

G(h,xi,yi; . . .;t N ,x N ,y N ) 

I) §{t x ,xx,yx) ■ ■ ■ ®(t N , x N , y N ) |0) 

via the LSZ reduction formula. The time evolution of the 
state is generated by the Hamiltonian 



// = / dx / dy ( in 2 + l -(d x <s>) 2 + i(a$) 2 + v@) 



where n = 9t$ is the conjugate momentum. 

Let us suppose now that the theory is restricted to the 
half space x < and the boundary condition is given by 
specifying some boundary potential at x = 

V B (<P(0,y,t)) 



The essential observation is that we can have two differ- 
ent Hamiltonian descriptions of this system. We can take 
t as the time variable and then the boundary is situated 
in space (i.e. it has a spacelike normal vector). In this 
case the Hilbert space consist of multi-particle states 



Hi 



= {\k±,%\\;k' x ,k'\\;...) B j 



where we indicate explicitly that the full momentum vec- 
tor is composed of two components, one perpendicular 
and the other parallel to the boundary. The states are 
normalized in the following way: 

B (k_L> fc|||fc-Li k\\) l B — 

(2ir) D u(k x , k\\)8 (All - fcj.) S [D ~ 1] (*|| - 

The subscript B indicates these states satisfy the bound- 
ary condition. In the asymptotic past all the particles 
move towards the boundary: fc^ > 0. For t — > — oo they 
are separated far from each other and from the boundary 
forming the in state, which is a free state. For £ — > oo all 
the scatterings and reflections are terminated they move 
backwards from the boundary are far from each other and 
from the boundary forming the free out state. These two 
Hilbert spaces are connected by the reflection matrix 



°B ut H 



The simplest case describes the one particle elastic reflec- 
tion on the boundary and can be written as 

R{0 imcS {k\\)) {2n) D 8 {0 - 6') 5 {D -^ -fcf|) = 
k \ , fen 



out 
B 



(k ± , fc|||KLL> k\\!b 



where we exploited the unbroken Poincare symmetry (in 
the coordinates t, y) and parameterized the perpendicu- 
lar momentum as 



k± = m c ff(fc||) sinh6> ; m e ff(fc||) = 
and then 

5(0-0') = m eS (k\\) coshO 5 (k± - k' ± ) 

The reflection factors can be expressed in terms of the 
correlators 

G(h,xi,yi; . . .;t Nl x N ,y N ) = 

B (0| $(ti,xi,sfi) . . .<5>(t N ,x N ,y N ) \0) B 

via the boundary reduction formula Q . The time evolu- 
tion of the state is generated by the Hamiltonian 



H B = 



dx I dy 



-V($) + 5(x)V B ($) 
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In the second description of the boundary theory we can 
take r = ix as time and £ = —it as space variable in 
the Hamiltonian formalism The Hilbert space now 
is that of the bulk theory (without any boundary), and 
time evolution is given by the bulk Hamiltonian (this is 
the extension of the crossing symmetry of the bulk the- 
ory, and is analogous to the open-closed duality in string 
theory,). In this case the boundary condition appears in 
time and serves as an initial state in calculating correla- 
tors: 

G(n,£i,yi;...;TN,£ N ,y N ) = 

(0|$(7i,fi,ift)...*(7>r,ew,iw)|B) 

where now \B) is a state in the bulk Hilbert space 7i. 
The boundary state \B) is in fact defined by the equality 
of correlation functions calculated in the two pictures. 
Using asymptotic completeness it can be expanded in 
the basis of the asymptotic in states of the bulk theory 



\B) 



K 



d^k 
(2^ 



...}|0) 



where A'(0,k\\) is the in asymptotic creation operator 
and the ellipses denote the terms with higher particle 
number. Due to translational invariance in the spatial 
direction all the contributing states must have zero to- 
tal momentum. Ghoshal and Zamolodchikov have shown 
that the coefficient K 2 is related to the one-particle 
elastic reflection factor R on the boundary by 



K 



(Q,m cS (k\\) 



R 



- 9,m eS (kn) 



(here we simply treat ku as a label for infinitely many 
species of two-dimensional particles of mass m e s(k\\)> 
using the fact that fci is conserved in the reflection off 
the boundary). The one-particle coefficient K\ is only 
nonzero when the vacuum expectation value of the inter- 
polating field $ of the particle does not vanish. In this 
work we suppose that b (0| <& (t, x, y) |0) s = (which is 
the case for the usual applications). We plan to present 
the general case in a forthcoming publication, together 
with more detailed derivation of the form of the coeffi- 
cient Ki and K 2 using the boundary reduction formula 
obtained in ;6], and a generalization to more complex 
particle spectra with several different masses. Here we 
restrict ourselves to a vanishing K\ to keep the discus- 
sion simple. We remark that when the theory in the bulk 
is free and the reflection is elastic, the boundary state can 
be written in a closed form 01 



\B) 



exp 



dO 
2^ 



iD-U 



,D-1 



(27T) 

K 2 (9, m off (fc| 1 )) Al(-9, -k 1 )A\ n (9, k 1 ) 



(2) 



We note also that if we have more then one particle 
species, with the same mass m created by A' in (6, h\)j 
(i.e. a multiplet), the formula for the boundary state 
changes as 



\B) 



00 d9 
2^ 



(27T) 



D-l 



K$ (9, m eS (A, 1 )) At (-9, -k\ 1 )iA\ n (9, fc, , ) + 



|0) 



Boundary conditions considered in the context of the 
Casimir effect generally allow transmission as well, and 
such boundaries are called 'defects'. A suitable gener- 
alization of the above formalism can be obtained by a 
folding trick, which maps the defect into a boundary sys- 
tem jjj. Suppose now that a defect is located at xq. In 
the crossed channel picture it can be represented by a 
defect operator which acts from the bulk Hilbert space 
of the x < xq system into that of the x > xq system. Lets 
denote the operator creating the particle for the x < Xq 
domain as A\ while for the x > xq domain as A\. There 
are now four one-particle reflection amplitudes: R are 
the ones preserving the species number 1,2, while 
are the ones changing 1 into 2 and 2 into 1, respectively. 
Upon the folding correspondence, R correspond to re- 
flection processes on the two sides of the defect, while T ± 
describe the transmission amplitudes from one side to the 
other. Using the folding map to the boundary system we 



ip 1 

obtain the defect operator [7| as 



D 



1 + 


/•OO 

J 00 


R+ 


/iir 

VT 


T+ 


/in 

VT ' 


T~ 


/Z7T 

VT 


R- 


/in 



d9 f dP-ije 



4vr7 (2tt) 



D-l 



(3) 



T+[— - 9,m cS (k ll )\A\(-9,-k ll )A 2 (-9,-k ll ) + 



m^mMAK-e^k^AKe,^) 



-e,m e s(k l] ))A 2 (e,k {l )A 2 (-e,-k ll ) + 



which (with the same conditions as for the boundary 
state, but now elasticity is required for the combined one- 
particle reflection/transmission amplitude [8j |19|) can be 
exponentiated to the form 



D 



exp 



dB_ r d D - 1 k 
4vri (2tt) 



D-l 



(4) 



R+(^-e,m cS (k ll )y\(-9 : -k n )A\(9,k ll )+ 
T+(y - 0,m e ff(fc||))4(-0,-£||)^2(-0,-fc||) + 



,m eB {k {{ ))A 1 (e,k l{ )A^{0,k { \) 



. 1 



+R- - 9, m eff (fc||)J A 2 (9, k l{ )A 2 (-9, -fc N ) 
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B. Derivation of Casimir energy 

Let us now turn to the derivation of Casimir energy of 
a D + 1 dimensional scalar field <&(t, x, y) in a domain of 
width L in x. To facilitate later applications, it is useful 
to consider the case when the field $>(t,x,y) is allowed 
to penetrate through the ends of the domain in such a 
way, that in the two domains of width R adjoining L, 
the dispersion relation of the asymptotic particles may 
be different from the vacuum one in L. The totally re- 
flecting boundaries at x = ±(R + L/2) are treated as 
auxiliary boundary conditions which are necessary in or- 
der to have a discrete spectrum; at the end of calculations 
we shall take the limit R — > oo and check that the results 
are independent of the choice of the auxiliary boundary 
conditions. So let us take the following situation in the 
coordinate x 




where total reflection occurs at xi = —R — L/2 and 
x r = R + L/2 with reflection factors Ri(9, m e ff(fc|i)) and 

Rr(6,m e ff(k\\)) and there are two "defects" (i.e. bound- 
aries allowing transmission) located at x\ = —L/2 and 
X2 — L/2 with one-particle defect matrices 



Di(9,m eS (k\\)) 



R+(0,m eS (k {l )) 
T+(0,m eff (£||)) 



"(0,m e ff(fe||)) 
~(0,m e ff(fc||)) 



where i = 1,2. The Hamiltonian of the system is denoted 
by Hb and its Hilbert space by Tt b ■ It is the ground state 
energy of the system which is of interest for calculating 
the Casimir effect, i.e. the lowest eigenvalue of Hb, which 
can be evaluated via the partition function. Compacti- 
fying all infinite (temporal and spatial) dimensions on 
circles with perimeter T we can calculate the partition 
function of our system in two different ways: 



Z R (L,T) = Tr HB e- THs 

-B.Hi 1 '' 



(Bile" 



Z?ie 



\B r ) 



where is the periodic Hamiltonian in the x channel 
in the three domains indexed by i — 1,2,3. Inserting 
complete sets of bulk states, taking the limit R — ► oo and 
normalizing the ground state energy in infinite space to 
Eq = we obtain 



Zoc(L,r) = ^2{B l \0){0\D 1 \n){n\D 2 \0){0\B r )e 

n 

= £<0|I>i|n><n|I> 2 |0>e 



— LE n 



-LE„ 



where the dependence on the auxiliary boundary condi- 
tions Bi r drops out (since in the R — > oo limit only the 
vacuum state contributes from the expansions of \B^ r )). 
The first few terms can be written explicitly as 

l + EE (0|£>i|^fcn;^,?i|)(^fc||;^gl||£>2|0) x 



e,k\\ e 'A. 

e -L(m otf (fc||)coshe+m cff (^|)cosh0') _|_ q/ 



-3mLj 



The term 1 is the contribution from the vacuum (|n) = 
|0)), the next term comes from two-particle terms in (|3J) 
and the higher-order corrections come from the higher 
multi-particle terms. This is a sort of cluster expansion 
similar to the one used in Q, valid for large values of 
the volume L. It is obvious from these expressions that 
the leading (two-particle) contribution depends only on 
R^ and i?^. The ground state (Casimir) energy (per 
unit transverse area) can be extracted from the partition 
function as 



E(L) 



lim 



>(L,T) 



The result is 



E(L) = 



(16 

47T 



cosh( 



R^\^— + e,m eff (k ll ))Rj[— -0,m eff (fc||)J x 

e -2m off (fc||)cosheL _|_ 0( e ~ 3mL ^ (5) 

The correction terms correspond to higher particle terms 
in the expansion Q of the defect operator D and in- 
clude the amplitudes of reflection/transmission processes 
involving more than one particle in at least one of the 
asymptotic states. These can be computed e.g. using 
a BQFT formulation as the one in 0, but it is obvious 
that they are suppressed by a factor e~ mL with respect 
to the leading order term due to the presence of at least 
one additional particle in the corresponding term of the 
expansion of the defect operator D. 

The formula JSJ, describing Casimir effect to leading 
order at long distances, is the main result of this paper. 
Note that it is applicable in the presence of nontrivial 
bulk and boundary interactions: their effects at leading 
order are contained in the reflection factors R^ , so as long 
as there is some theoretical or experimental input from 
which these can be determined the leading order contri- 
bution can be evaluated. In integrable (2D) boundary 
theories are obtained as solutions of the boundary 
Yang Baxter equation, and the bulk interaction manifests 
itself thorough the bulk S matrix appearing in this equa- 
tion. In case of nonintegrable interacting bulk theories 
e.g. a perturbative expansion can be given to determine 
R ± , (for more details see jy|). 

Another important point is that this approach formu- 
lates the Casimir effect from an infrared viewpoint. Stan- 
dard derivations of the Casimir effect (such as the one 
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presented in the appendix) solve the microscopic field 
theory. This necessitates tackling diverse issues such 
as renormalization, and also the possibility that the in- 
frared (long distance behaviour) may be quite different 
from the microscopic description of the theory (as is the 
case for example in QCD). Formula is expressed in 
terms of the asymptotic particles, the long distance de- 
grees of freedom, and provides a long distance expansion 
for Casimir energy. It describes a finite size effect in a 
boundary quantum field theory close in spirit to our pre- 
vious investigation of the boundary Luscher formula , 
and it is not difficult to see that the latter is just a spe- 
cial case of |[EJ . As a consequence our formula has already 
been tested in interacting (integrable) 2D quantum field 
theories in (This is the only case to our knowledge 
where exact reflection factors of interacting BQFTs have 
been computed). Furthermore, the result (jSJ) includes 
the contribution of states localized to the defects (called 
'surface plasmons' in the literature), since they are taken 
into account as poles at imaginary rapidity of the reflec- 
tion factors. This fact is also demonstrated explicitely 
using the zero mode summation method in Appendix A. 

A very appealing property of the boundary state ap- 
proach is its universality and we shall see in the next 
section that it indeed reproduces all the results previ- 
ously known for the planar situation. But in those cases 
we can go even further, because in most calculations of 
the Casimir effect the bulk is free and the boundary scat- 
tering is elastic, therefore one can use the exponentiated 
boundary state @ together with the analogous defect 
operator Q which makes it possible to sum up all the 
multi-particle terms of the cluster expansion. The calcu- 
lation to be done is essentially identical to the derivation 
of the boundary Thermodynamic Bethe Ansatz equation 
with trivial bulk S matrix and reflection factors and 
R 2 , which was already performed in 0- The only dif- 
ference is that supposes fermionic statistics and that 
in our case the reflection factors R± and R 2 are not uni- 
tary due to the existence of transmission, but these do 
not change the overall reasoning of the derivation. The 
result is 



E(L) = ± I — cosh 

. 4-7T 



(^T™e ff (fc||) 

/ / -7T _f \ / —, 

l0g( 1 T R l[-J + 0."lcff(fc||)J.ffij"(y - 0,TO e ff(fc||) 



-2m c ff | ) cosh OL 



(6) 



(the upper/lower signs correspond to the 
bosonic/fermionic case, respectively). Although formula 
is more general and adequate enough for the large 
distance regime, the theoretical situations considered 
in the literature can be directly compared to ©. It 
is obvious that for L mT 1 © reduces to (JSJ. The 
Casimir force is dominated by modes with momentum of 
the order 1/L ; for large enough separation the bulk in- 
teraction can be dropped due to its short-ranged nature, 



and at the same time the characteristic energy is lower 
than the threshold for inelastic boundary scattering. 



III. APPLICATIONS 

To show the universal nature of these results we con- 
sider some applications. We demonstrate that taking the 
reflection factors from the literature and substituting into 
immediately gives the Casimir energy without going 
through a detailed analysis of the microscopic quantum 
fluctuations. 

As a first application we determine the Casimir energy 
for a massive free scalar field, V(Q) = n ^-^> 2 , subject 
to Robin boundary conditions, Vb(&) = f^ 2 , on two 
parallel hyperplanes. If the two planes are located at 
x = and x = L the boundary conditions are given as 

<9 a $ - C!$| x= o = 0; 9 X $ + c 2 $U=z, = 0; c u c 2 >0, 

and the reflection amplitudes on these planes can be writ- 
ten as 



Rj(9,m eS (k\\)) 



m e ff(A:||) sinh# — icj 



J = 1,2. 



fTfieS (fc||) sinhfl + icj 

Note that they have the same form as in the two di- 
mensional case, but now they depend also on fen via the 
rapidity parameterization k± — m c ff (fc||) sinh(9. When 9 



is continued to 9 — > 9 + i~, we have 



+ i^-,m eS (kn)) = 



m c fj(fc||) cosh# — Cj 
m e ff (fc| | ) cosh 9 + cj 



.7 = 1,2. 



Introducing the variable q by m c g (fc||) cosh# = 
y/ m 2 + q 2 and performing the angular integrations one 
obtains 



E(L) 



(Att) d / 2 T(D/2) 



dqq 



D-l 



log ( i - ^ m2 + g2 ~ Cl vjS + q2 ~ C2 C -2^A^ I _ 

\ yjm 2 + q 2 + ci m 2 + q 2 + c 2 J 

This formula is a new result of the present paper. To 
show its correctness we consider two already known lim- 
iting cases. First consider the limit when both Cj — > or 
Cj — > oo corresponding to Neumann or Dirichlet bound- 
ary conditions for In both cases the coefficient of 
the exponent in the logarithm becomes one, and we ob- 
tain the Amb jorn- Wolfram result for Dirichlet boundary 
conditions as reported in Milton's book 0- In the 
second limit we let the mass of the scalar field vanish 
to — > 0; then 



E(L) 



(4tt) d / 2 T(D/2) 



d g ^log(l- g - Clg - C2 e-^ 
q + Ci q + c 2 j 
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It is straightforward to show that this result coincides 
with that obtained in 11] for the Casimir energy per 
unit area of a massless scalar field with Robin boundary 
conditions. 

The second application concerns the Casimir force be- 
tween parallel dielectric slabs separated by a vacuum slot 
of width L, a problem first investi gate d by Lifshitz and 
collaborators almost fifty years ago [lj ■ In this geometry 
the dielectric constants in the three regions are: 

e(x) = e%; x < 0; e(x) = £2; L < x; 

e(x) = 1; < x < L. 

In this case Rf , 2 (i^ +$) are nothing else but the (appro- 
priate analytic continuations of the) ordinary reflection 
amplitudes of electromagnetic waves incident from vac- 
uum at the plane interface between the vacuum and the 
dielectric materials, given in many textbooks (e.g. [l3|). 
Indeed denoting again by (k±,ku), 10 the wave vector and 

frequency of the electromagnetic radiation (fc|| being the 
component parallel, while k± the one perpendicular to 
the plane interface) the reflection amplitudes are 



; 2 - kl - Je,cu 2 - kl 



; 2 - kf, + \ UiUj 2 - kl 



^perp^fcll) 



when the electric field E is perpendicular to the plane of 
incidence, and 



11 



£i(w), 



when E is parallel to the plane of incidence and i = 1,2. 
(To obtain these expressions we assumed that the per- 
meabilities of the slabs are unity, fii — 1). Continuing 
8 — > 9 4- corresponds to continuing to to purely imagi- 
nary values oj — > 77i ff(fc||)isinh# = iQ, and the reflection 
amplitudes become (q 2 = kf<): 



V( 2 + v 2 + v^C 2 + q 2 



- ^/aC 2 + q 2 



iVe + ^+V^C + q 2 



Since the electromagnetic field can have both polariza- 
tions, one obtains the following form for the Casimir en- 
ergy per unit area 



1 



E(L) = — I d q z I d( 



log (l - flWpCC, q)R^ p (i£ q)c- 2L V^ 
log (l - i?« (iC, q)R$M, q )c- 2L ^^) 



The Casimir force ^{L) = —dE(L)/dL computed from 
this expression agrees with that of Lifshitz et al. as re- 
ported in the reviews 0, Q ■ 

Next we consider the Casimir energy of a massless 
fcrmion field in 1 + 3 dimensions subject to the "bag 
boundary condition" in our planar geometry: 



(1 - i-f 3 )i>\ x =o = 



(i + i7 3 )^U=L = 



rality but commute with E 3 = 



Furthermore 



Using the chiral representation of Dirac matrices one 
readily shows that these boundary conditions break chi- 

<j 3 

CT 3 

they imply that for both the positive (u) and the nega- 
tive (v) eigenvectors of £3 the reflection amplitudes are 
constants, satisfying R~^(in/2 + 9)R 2 ^(iTr/2 + 8) = — 1 

for s — u,v. Thus (taking into account crossing-unitarity 

00 

E(L) = -2^L J dqq 2 \og(l + e- 2L "), 


which indeed coincides with the known result 0, Il4| . 

IV. CONCLUDING REMARKS 

The main result of the paper is formula (0 which ex- 
presses the Casimir energy in large volume in a general 
form valid for any QFT - even for theories interacting in 
the bulk - in terms of the reflection amplitudes (infrared 
data). It passed the test of interacting but integrable two 
dimensional quantum field theories in Q . The boundary 
state approach provides a systematic infrared (large vol- 
ume) expansion, which in the case of noninteracting bulk 
and elastic boundary scattering can be summed up lead- 
ing to ©• In order to check the result for noninteracting 
bulk we present in the appendix an alternative derivation 
using the naive mode summation method. 

In the usual calculations of the Casimir effect correc- 
tions for temperature and roughness effects must be made 
before comparing to the experiment. Both corrections 
can be introduced into lloKit in a straightforward man- 
ner, followin g th e usual procedure in the literature (cf. 
@) see also |l6j . In the case of a massive theory with 
mass gap m the terms coming from higher particle cor- 
rections to £[J and © are suppressed by exp(— mL). For 
theories with a zero mass gap the situation is different 
and indeed the asymptotic state formalism is not well- 
founded if interaction with massless particles (such as 
photons) is taken into account, due to infrared divergen- 
cies. In the case of electromagnetic field the corrections 
to JSJ arc known to be suppressed by a\/L where A is 
the Compton wave length of the electron and a is the 
fine structure constant (Chapter 13 of Q), due to the 
fact that the self-interaction of the electromagnetic field 
arises only through radiative corrections and the leading 
one is an electron loop. 



7 



The striking fact about the general formulae l">l(ili is 
that they are manifestly finite, and only depend on phys- 
ical quantities, such as the dispersion relation of asymp- 
totic particles and their reflection amplitudes on the 
boundaries. This is even more emphasized in the bound- 
ary state formalism, where the bulk and boundary di- 
vergences appearing in usual derivations of the Casimir 
effect are manifestly absent, and the whole derivation is 
performed in terms of finite (renormalized) physical ob- 
jects. The appearance of divergences in the naive method 
of summing zero-point energies is a general feature when 
one calculates the interaction energy of point-like sources 
by integrating the energy stored in the field generated 
by them. The divergences are then eliminated by notic- 
ing that they contribute to the self-energy of the sources 
themselves and are present even at infinite separation. 
These divergences are absorbed by renormalization, but 
they are entirely unphysical and eliminated when ex- 
pressing the measurable energy differences in terms of 
physically meaningful, finite quantities. Further bonus 
of the general formulae (loKiH is that the reflection factors 
appearing in them can be calculated easily in the semi- 
infinite geometry, e.g. using perturbation theory along 
the lines of [||. 
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Appendix A: CASIMIR ENERGY BY MODE 
SUMMATION 



Here we want to support our derivation coming from 
the boundary state formalism by obtaining the analo- 
gous result using a new version of the mode summa- 
tion method for fluctuating fields which are free in the 
bulk. Ideas somewhat similar to these have been used 
in a similar framework earlier in |15| and their find- 
ings also support our general expression. The novel 
viewpoints of our approach are twofold: first we con- 
centrate only on the finite size piece of the ground 
state energy that should be independent of regulariza- 
tion and renormalization thus we do not need to dis- 
cuss the regularization/renormalization of the (infinite 
space) bulk/boundary ground state energies, since even- 
tually these are subtracted. The second feature of our 
approach is that we do not solve the quantization con- 
dition for the zero point frequencies explicitly but use 
only the functional form of these equations to obtain a 
general expression for the ground state energy containing 
only the reflection amplitudes. 

For simplicity let us start with the 1 + 1 dimensional 



free scalar field x) of mass m in a domain described 
earlier. We assume that in all three domains the fre- 
quency of the modes of $ is the same, but the wave vec- 
tors k', k" may be different from k (e.g. the mass of the 
scalar field is different in the three domains), but they 
can be regarded as functions of k. All the wave number 
parameters k, k! ', k" are chosen positive. Then for the 
plane wave modes one has the situation displayed in the 
following figure 




where total reflection occurs at x\ = —R — L/2 and 
x r = R + L/2 with reflection factors Ri and R r 

Ae - ik >(R+±) = Be ik '^ R+ ^Ri(k') (Al) 
Fe -ik"{R+k) = Ee lk " {R+ ^R r {k"), 

and the two "defects" (i.e. boundaries allowing trans- 
mission) located at x\ = —L/2 and X2 — L/2 describe 
one-particle scattering with defect matrices 



Di(k) 



Rf(k) r-(k) 

T+(k) R~{k) 



1,2 



In the case of elastic defect scattering their is no particle 
production and the defect matrices are unitary. They 
connect the incoming and outgoing amplitudes as follows 



Ae -ik'L/2 



£, e +ikL/2 



Be +ik'L/2 

De ^kL/2 \ ^ / Ce+* kL ' 2 

Ee +ik"L/2 J - 2 '"7 I p c -ik"L/2 



(A2) 



Unitarity and time-reversal gives the relations 

A(*) f A(A!) = 1 , Di{k)- 1 = Dii-k) 
Rr,i(k)*Rr,i(k) = l , Rr^k)- 1 = R rA (-k) 



(A3) 



Consistency of the homogenous linear system (|A1IA2|) is 

expressed by the quantization condition 



Q(k n ) = 



(A4) 



where 
Q(k) = ( I 



Ri(k)Rt(k)e 



ilk' R 



R r {k)R2{k)e 



i2k" R 



e' 2kL R^(k)R+(k)[ 1 



Ri(k)e 



i2k' R 



R+(-k) 



1 - 



Rr(k) 



ilk" R S 



Rz(-k) 
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and we exploited l|A3(l to write 

det A - Rl{k)/Rj{-k) = R+(k)/RT{-k) 

Due to unitarity the quantization condition can only have 
real and purely imaginary solutions. For each real solu- 
tion k n , —k n is also a solution, and fco = always solves 
(|A4I) . The purely imaginary solutions located at kj = inj 
(0 < Kj < m) are related to the poles of the reflection 
and transmissions factors and correspond to defect and 
boundary bound states. Then the ground state energy of 
the system is given by the sum of zero modes, which can 
be written in an integral form 



E(L,R)= \"{kn) + Y,\^ iK i 



{fc„>0} 




1 \- 

2 1 2 



where Q'(k) — oj(k) = y/k 2 + m 2 , h is set to one and 
the contours are shown on the figure 



pole 

( 

pole 



iK 3 ) c 3 



iK,] C 1 



-k 3 -k 2 -k! o kj k 2 k 3 



The poles on the figure denote the poles of the Q func- 
tion originating from the poles of the reflection and trans- 
mission amplitudes. Contrary to Kj, which do depend on 
the volume, the location of these poles is independent of 
L, R. 

We concentrate on the L dependence so the L inde- 
pendent Co integral can be dropped. Next we analyze 
the contribution of C_ . Using i|A3(l we can relate 



n-k) Q'(fc) 



Q{-k) 



Q{k) 

l0gi?r 



= -2iL-2iR 



dk' 
~dk 



dk 



dk 



log(detL»i_D 2 



dk" 
~dk 
d 

~ dk 



log Ri. 



The first two terms on the right hand side (proportional 
to the volumes L, R) correspond to divergent bulk en- 
ergy contributions. The other terms here correspond 
to infinite boundary and defect energies which are in- 
dependent of the volumes (L, R). The precise defini- 
tion of these terms needs some regularization and renor- 
malization, but they give no contribution to the Casimir 
force since they are present in the infinite system as well. 



Therefore they can be dropped and then the C+ and C_ 
integral terms give identical results. 

Thus the relevant contribution to the Casimir energy 
can be written in terms of the variable m sinh 9 = k as 



E(L, R) 



00+111 d0 uaQ'W 
— cosh e — — 

-oo+i V 4th Q{6) 



(where 77 is a small positive parameter). Next we shift 
the contour to 77 = -| which eliminates the discrete sum 
via the contributions from the purely imaginary zeros of 
Q(0). There can be additional contributions from poles 
of Q (6) on the imaginary axis, but their positions are 
independent of the volume and therefore they only con- 
tribute volume independent terms which can be dropped. 
Integrating by parts and introducing the new variable 
9 -> + 9 we obtain 



E(L, R) 

In the limit R — > 00 
d9 

47T 



— cosh# logQ( — 
4tt WV 2 



E{L) 



cosh( 



log 1 - i?r 



9 Ri 



Y 



-2m cosh 9L 



It is straightforward to generalize these expressions to the 
D + l dimensional case when the boundary conditions are 
specified on two parallel D — 1 dimensional hyperplanes. 
Introducing the component parallel (kn) respectively per- 
pendicular (k±) to the hyperplanes one can write 



k = (fc_L,fc||); uP- - k\ = 
k± — m e s(kn) sinh^. 



kt 



'■err 



(fcn); 



(A5) 



We assume that the system is translationally invariant 
in the direction of the hyperplanes, and so the reflection 
and transmission processes preserve fen. Therefore we 
can apply the previous considerations to the decoupled 
one-dimensional systems labeled by the parameter k\ \ and 
sum up their contributions. Note that both u> and the 
reflection/transmission amplitudes depend on kn, but we 
shall not write this out explicitely. The Casimir energy 
per unit transverse area can be written as 



E(L,R) 



(A6) 



j-D-li 



— m off (fc||) 



^ cosh 6 lo § q(y 



From this, the Casimir energy between the two defects 
can be obtained by taking the limit R — > 00 



E(L) 



d 



D-li 



'II 



(2 7 r)^i Wcff(fc " ) l 00 ^ COSh 



d6 



log ( I - R- [ - 



(A7) 

e -2m ff(fc||)cosheL \ 
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Note that the latter result is independent of the auxiliary 
boundary conditions R r and Ri which can be expected 
on physical grounds. Taking into account the crossing- 
unitarity property satisfied by the reflection amplitudes 
0,0 which for free bosons/fermions takes the form 




one can see that this result is in a complete agreement 



with the formula obtained form the boundary state for- 
malism previously ©■ (The results in l|A6l IA7|) should 
be multiplied by (—1) in case of computing the Casimir 
energy of fcrmionic fields, since their vacuum energy is 

-££w(k)). 

If there are more than one fields in the problem, or the 
fields have more than one component one has to sum up 
the contribution of all the fields (components) to obtain 
the Casimir energy. 
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